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FOREWORD 


This report describes a portion of the results obtained on 
NASA Grant NSC 3044* TUIb work was done under subcontract to the 
University of Illinois, Urbana, with Prof, S.S. Wang as the Principal 
Investigator. The prime grantee was the Massachusetts Institute of 
Technology, with Prof. F.J. McGarry as the Principal Investigator 
and Dr. J.F, Mandell as a major participant. The NASA - LeRC Project 
Manager was Dr. C.G, Chamis. 

Efforts in this project are primarily directed towards the de- 
velopment for finite element analyses for the study of flaw growth and 
fracture of fiber composites. This report presents exact solutions for 
edge delaminations in angle-ply composites. 


ABSTRACT 


Edge delamlnatlon has caused severe concern in the design and 
analysis of advanced composite materials and structures. Due to its 
complex nature , very' limited knowledge for the problem is currently 
available. It involves not only geometric and material discontinuities 
but also inherently coupled mode I, II and III fracture in the layered 
anisotropic system. Based on comp lex- variable stress potentials in 
the anisotropic elasticity theory and eigenfunction expansion* exact 
orders of the crack-tip stress singularity and complete field solutions 
are obtained. Results are given for edge-delaminated composites 
subjected to uniform axial extension for illustrative purposes. 

Effects of geometric, lamination, and crack variables are determined. 
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1* INTRODUCTION 


' ) 

Edga delamination Is frequently encountered in angle-ply compo- 
site laminates. It is due to high stress concentrations at geometric 
boundaries and the inherently weak interlaminar strength along the ply 
Interface, The problem posed by edge delamination is of great theoreti- 
cal interest* It also is of significant technical importance in determining 
the structural integrity and damage tolerance of advanced fiber-rein- 
forced composites and their applications to advanced engineering 

structures and components. The presence and growth of delamination 

\\ 

cracks from geometric boundaries of composite laminates may lead to 
severe reliability and safety problems of fiber composite materials 
and structures such as the reduction of structural stiffness, the 
exposure of the interior to adverse environmental attack, and the 
disintegration of the material, which may cause the final failure, 

Thus, understanding the basic nature of edge delamination is of 
critical Importance in damage characterisation and accurate assessment 
of flaw criticality and structural integrity of advanced composites. 

The edge delamination problem is very complex in nature and 
extremely difficult to solve . It involves geometric and materials 
discontinuities, i.e. , free edges, interlaminar craqks and variation of 
ply properties in the transverse direction. It also Involves inher- 
ently coupled mode I, II and III fracture In the anisotropic layered 
material system such as angle-ply composite laminates. Edge de- 
lamlnatlon is basically a fracture problem involving an interfacial 
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crack between two highly anisotropic fiber-composite laminae 
under general loading conditions. The problem of an Interfaclal crack 
between two dissimilar isotropic materials has received much atten- 
tion recently, for example, Refs. (1-12], But the study of a de- 
laminatlon between strongly anisotropic fiber composite layers, 
especially In finite-dimensional laminates under general loading, 
has been very limited, to the author's knowledge. In this paper, a 
rigorous investigation of the basic nature of the coupled opening, 
sliding and antiplane shearing fracture behavior of edge delamination 
is presented for composite laminates under uniform axial extension. 
Basic formulation of the problem based on the theory of anisotropic 
elasticity and eigenfunction expansion is given in the next section. 
General solutions and associated stress singularities for the edge 
delamination are derived in Section 3. Fracture parameters such as 
mixed-mode stress intensity factors Kj., and K IIX as well as the 
energy release rate G are defined and determined in Section 4. Num- 
erical results for edge delaminated composite laminates subjected 

to uniform axial loading are shown in Section 5 to illustrate the 

■j 

fundamental behavior of edge delamination cracks. Effects of geo- 
metric, lamination and crack variables are studied in detail. 


d If 
ij 
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2. FORMULATION 


Formulation of the delamination problem is baaed on the theory 
of anisotropic elasticity for nonhomogeneous solids* Well known 
Lekhnitskii'a stress functions {133 are introduced to establish 
governing partial differential equations for field variables. An 
eigenfunction expansion method is employed for determining the stress 
singularity at the delamination crack tip. The boundary collocation 
technique is then used to evaluate the complete solution for finite- 
dimensional composite laminates, 

2.1 Assumptions and Delamination Model 

Consider a composite laminate (Fig, 1) composed of unidirec- 
tional fiber-rcinforced plies of uniform thicknesses, h^, h 2 ,..,.h n . 

» 

For simplicity but without loss of generality, we restrict ourself to 
the cases of symmetric, angle-ply composite laminates with fiber 
orientations (O^/O^/* • • ^cknesses are also assumed to 

be symmetric with respect to the x-z plane, i.e., h^ * h , h 2 * h n _^, 
.... The composite has a finite width 2b and is subjected to a uniform 
axial extension, e, where e - constant, along the z-axis. The composite 
laminate is sufficiently long that, in the region far away from the 
ends, end effects are negligible by virtue of the Saint Venant 
principle. Consequently, stresses in the laminate are Independent of 
the z-axis. The case in which stresses and displacements are independ- 
ent of z corresponds to the well known generalized plane deformation 
[ 13 ] , Edge delamination occurs in the form of a crack along the 
interface of dissimilar plies with fiber orientations, 0 k and 9 k+1 . 
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Perfect bonding is eesuned in the composite everywhere except the 
region of delaainetion. 

2.2 Basic Equations 

For eech individual lamina* the constitutive equations in the 
structural coordinates x-y-z may be expressed by the generalized Hooke's 
law in the contracted notation as 


c i - s u °i 


(i*j ■ 1*2*.. ..6) 


( 1 ) 


where and are strain and stress tensors* and S^* the compliance 
tensors, respectively. The strain-displacement relationships are given 
by 

(2a-b) 


e l " e x " u ’x> 


e 2 " e y " v *y’ 


(2c-d) 
(2e-f) 

where the subscript after comma denotes the partial differentiation with 
respect to the variable. Under these assumptions, there remain three 


■ e ■ w, , 

e. « y 
4 T yz 

- w*,. 

+ V, » 

z z 

y 

z 

■ Y ■ W, + U, , 

’xz ’x *z’ 

c- ■ Y 
6 'xy 

m u, 

+ v, * 

y 

’x 


compatability relations: 


if 


+ E y, ■ Y xy, ’ 


’yy Y ’xx ’xy 

^" Y xz, y + Y yz, x ) *y " ° 9 

(“Y V2 _ + Y*, >»„ - 0. 


'yz, x ’ 'xz fy ’x 


(3a) 

(3b) 

(3c) 


For a symmetrical angle-ply composite laminate subjected to uniform axial 
strain a, it can be shown easily that 
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0, 


(4) 


» 


Y - Y 
XZ.y yz. x 


sine® the relative angle of rotation for the symmetric composite lamin- 
ate about s«axis vanishes. 

Based on the definition of th<i problem, i.e., e z - * - constant, 
o z in Eq 1 can be expressed in terms of other stress components by? 

a z m " s 3j a j^ s 33 (J " 1)2, 4, 5, 6). (5) 


Thus, the generalized Hooke’s law may be modified to have the following 
form: 

e i " ^ij °j + e l 

where 

®ij " S ij " S 3lSj^ S 33’ e i " eS i3 /S 33 ** 3) ’ 

It can be seen from Eq 6 that e i has the role of initial strains in the 
laminate, For the current stress formulation, it may be more convenient 
to Introduce the initial stress a^ Q such that 

e i * ( a j " a jo )f (i,j * 1,2, 4, 5, 6), (7) 

~-l 

where o iQ ■ -S^ e^ , It is possible to decompose the complete solutions 
into two parts, i.e., 

(8a) 
(8b) 

e r > * *ij Co i P) - V ’ (9a ' b) 


w _ ^(h) + (p) 

a i " a i C J 


e (h) + c (p) 
e i * e j 


where 


- S,. and 

i 1 j J 


(6a) 


(6b) 
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3, GENERAL SOLUTIONS AND ASSOCIATED STRESS SIMILARITIES FOR DELAMINATION 


3*1 Solutions for and 

Introducing the. stress functions F(x,y) and T(x,y) which satisfy 
equations of equilibrium identically and following the procedure by 
Lekhnitskil [13], we obtain a pair of coupled partial differential 
equations as follows: 

; L 4 F + L 3 Y « 0, (10a) 

* L 3 F + L 2 ¥ « 0, (10b) 


where 1 2> L^ and L^ are llnfei# differential operators of the second, 
third and fourth order, respectively, defined by 

g2 » g2 g2 

L 2 * S 44 ax 7 " 2S 45 3x3y + S 55 3? r * 

«** ft 3 ** ** ft 3 •» « ft 3 

l 3 ■ _S 24 5P - + <S 25 + S 46> 3777 ‘ (S 14 + S 56> W 

~ 3 3 

+ S 15 W ’ 

L 4 ‘ ®22 3X 5- ” 2 *26 sJsy + (2 ®12 + ti&iy 1 

a** . - 

“ 2S 16 3x3y T + S 11 a? 5 " ’ 


(10c) 


aod) 


(10e) 


Lekhnitskil [13] has shown that the general solution for Eqs lOa-lOb may 
be expressed as 


F(z k> " .1 W - nz w > ’ VW 


k-1 


(lla-b) 
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where Z^ - x + y^y{ the prime (*) denotes differentiation of the function 
with respect to its argument? and are the roots of the algebraic 
characteristic equation 


and 


with 


«' 4 (h)^ 2 ^> " * 3 0 * ) * 0 , 

n k ‘ 'W'W “ "W-W 1 

V> - S 55 u 2 - 2S„|i + S w , 

t 3 (n) * S 15 U 3 “ (S w + Sj 6 )v 2 + (S 25 + S^ 6 )u - § M , 

V 5 - S U M 4 - 2S 16 H S + (2S 12 + S 66 ))l 2 - 2S 26 W + S^. 

Introducing the following form for the function F^Z^) 

5+2 

W - C k Z k /tC5 + 2>C« 4* 1) ], 


(12a) 

(12b) 

(13a) 

(13b) 

(13c) 


(14) 


where and 5 are arbitrary complex constants to be determined later, 
we can obtain the stress and displacement expressions as follows; 


°k h) ‘ J x [ Vk Z k + c k+3 4 t) 

oi h) - l Wi + C V 4., t) , 


1 k k ’ k+3 k j 


k-1 


< l z -:E >Wk + c k + 3\^- 


yz 


k-1 


^ - J^VkVk + °k+3 ‘‘k \ ^k 1 * 


(15a) 


(15b) 


(15c) 


(15d) 
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-(h) 

xy 



•""'Q 1 4 

. - u. Z .** 
k+3 w k k 


and 

“ <h> " Ji tCkP ^ +1 + °k+3 *k *k +1)/({ + 15 

VCH) “ Jx'WS* + C k+3 \ 5 t +1 ' /(S + « 

“ <h> " J 1 t Vk Z k +1 + C k+3 *k *k +lj /<4 + 1J 

where 

p k " ®ll M k * ®12 “ ^I4 n k + ®15 n k W k " ®16 y k 
q k " ®12 w k + ®22^k " ®24 n k^k + ®2S\ ” ^26 ' 

\ “ ®l4 M k + ^24 /M k " ^44 n k /M k + ®45 n k " ®46 * 


(13t) 

(16a) 

(16b) 

(16c) 

(16d) 

(I6e) 

(16f) 


The constant, 6, in Eqs 15 and 16 may be chosen that the stresses 
and displacements and u^*) satisfy interface continuity and 

homogeneous boundary conditions. Taking complex conjugate of Eqs 15 
and 16, their forms are invariant; thus, 6 appears as a set of complex 
conjugates, which enables to make Eqs 15 and 16 real functions by super- 
position. Furthermore, finiteness of displacements at the origin re- 
requires that Re[<$) > -1, where Re represents the real part of 6 . 

3.2 Solutions for <j£ p ^ and u^ 

Since cr io and e.^ in Eqs 6 and 7 are constant, we may choose er| p ^ 
in Eq 9 as constants so that they satisfy the equations of equilibrium 
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and compatibility conditions identically. Tor each individual lamina, 
let take the following form: 


4 P> - V + + \ ^ - 


(p) 


+ l Wk + V» 


yo 


k-l 


(17a) 

(17b) 


■ Vo ’ J^Vk + ^V* 


(17c) 


T xz' - T x» + ,.L (d kVk + \ \ V- 


(17d) 


k-l 


“ Vo" J 1 <Vk + d k ir k > - 


(17a) 


Substiutlng Eqs 17a-17e into Eq 9b and integrating the strain-displace- 
ment relations, we obtain 
3 


U<P> " , l, Wk + \ \ V - V + V + \ 

K" i 


(18a) 


V<P) ' Jl (d ^k Z k + d k “k Z k> + V ' V + V o' 


(18b) 


w 


(p) 


t < d k c k z k + d k *k V - V + V + + w c 

k-l 


(18c) 


where u , v , w and w, are related to rigid-body displacements and rota- 
o O O 1 

tions* The complex constants d^ are required to satisfy the near-field 
traction boundary conditions and continuity conditions along the ply 
interface. 
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3,3 De lamination Crack-Tin Stress Singularity 


Consider a delamination between two plies, say the kth and (k+l)th 
ply, in a composite subjected to general loading as shown in Fig, 2* 
Assuming that interlaminar crack surfaces are free from traction, we 
introduce the following boundary conditions for the eigen stresses 

.(h). 

tf i * 

erf 1 * *- T S‘ > " -0 (i - k, <t* « jr; 1 - k+1, * - -it). (19) 

y yz *y 

The superscript h is dropped in the expression for convenience. 

Continuity conditions for displacements and interlaminar stresses 
along the Interface, $ - 0* 


,„0O -<k) _(k), fo (k+l) (k+1) (k+1), 

' y ’ yz ’ xy 1 l °y ,T yz ,T xy ' 

(20a) 

<«< k \v <k >,.r (k >> - (.WD-^Wl^Wl),, 

(20b) 


Substituting Eqs 15 and 16 into 19 and 20, we obtain the following 

(k) (k+1) 

twelve linear algebraic equations in C^ and C^ 9 i 


l + o' 11 " 5 ^ r£>> ■ 0. 


m-1 


m+3 jm 


(J-1,2,3) 


(21a) 


l { Q “ ilT<s c ( k+1 ) r ( k+1 ) + T. (k+1) } - 0 (j-1,2,3) » (21b) 

m im m+o jm 


m-1 


l <*«? + w e x) > 


m-l 


m-1 


J m+3 rm 
(r-1,2,3,4,5,6) , 

pCk). r tk) n ( k ) P 0t>. u ( k > r (k). n Ck> r ( k >„ a C k ) and 

where r lm -1, - n m > ■ W m » r i™ " P m » Um and 


(21c) 


rP^- t^. Solving Eqs 21c for cf k ^ and substituting the resulting 
om m m 
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expressions Into Eqs 21a , we get 


l (C< k+1> l [e 1 ’^ r (k > - 

u n mn jm 


+ e 


n-1 


m-1 


a Y 1 > 

a (nrf3)n l jm 


0 . 


(21d) 


Equations 21b and 21d consist of a system of six linear homogeneous 
algebraic equations. For nontrival solutions of the determinant 

of coefficients of the algebraic equations must vanish. This leads to a 
characteristic equation of the following form: 

< e 12l,S - Dalits)! - 0, (22) 

where |a(<$)| is a 3 by 3 determinant involving 6 in a transcendental form 

and material constants, p ^ and p 'k*-) ^ (k+1) adjacent 

m m m m 

layers. Details of A (5) may be found in Ref. [14], The general form of 
6, which are the eigenvalues of the problem, may be written as 

6 n - n, or * (n - “■) ± iy (n-0,1,2, . . .) , (23) 

where y is a constant related to elastic constants of adjacent plies. 
Thus, for each 6^ we have the eigenfunctions of the form Eqs 15 and 16 
whose coefficients may be determined from the remote boundary conditions 
other than Eq 19. It is important to note that the <S n bounded by 

0 > Re[5 n ] > -1 (24) 

characterize the inherent stress singularities of the delamination crack 
stresses in a composite laminate. 

For cross-ply composite laminates, the differential operator 
vanishes identically. Thus, F(Z) and Y(Z) are uncoupled, and the form 
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of Eq 23 can be simplified and expressed explicitly as 


n 


n. 


or 


(n - |) ± An { [b + <b 2 


- 4a 2 ) T ]/(2a)}, 


(25) 


where a and b are related to material constants 

U< k > and shown in Appendix 1. In a limiting case of 

isotropic materials, it can be shown that -6 have the form. 


n 


Cn - |j ± Htl{[G (k) + G (k+1) (3 


4v <k) ) ] 


/ [G ( k+1 > + G (k) (3 - 4u ( ' k+1 h]). 


(26) 


where 0 and v denote the shear modulus and Poisson’s ratio, respectively. 
The eigenvalues of Eq 26 were first obtained by William [1] and later by 
Zak, at al. [2] for interfacial cracks In Isotropic media. 


4. DELAMINATION STRESS; INTENSITY FACTORS AND ENERGY RELEASE RATES 


The eigenfunctions and the unknown constants for and a^P) i n 
Eqs 15 and 17 are determined by imposing appropriate (materials and 
geometric) symmetry and traction boundary conditions, which will be 
discussed later. Hence, complete stresses and displacements °< o) and 
u<«> in the a-th lamina can be fully established. Neglecting the 
higher-order terms, we note that the typical structure of near-field 
crack tip stresses can be shown to have the following form: 


.. <*> - y y f £ (cO z 6 j + tt («) 
i ~ k i 1 [f ijk Z k + 8 ljk Z k ] » 


(27) 


where fj^O and g(<0 are re i atec } to the material constants, geometry, 

and boundary conditions; are eigenvalues bounded by 

-1 < Re[<S,] «i 0 to insure the positive definiteness of strain energy 

J 

of the elastic body. It is clear that the eigenvalues which satisfy 
Eq 24 lead to asymptotic near field stresses. For the convenience of 
further development, the stress a,£ a ^ around the crack tip may be re- 
written as 


" 1 ^fa^(x,y;6 , ) + 0 (non-singular, higher-order terms), 
1 j-i 


'1 


(28) 


(«) 


(a) 


where or£j 7 is the j-th singular component of the stress cor- 

responding to the eigenvalue <5^ which meets Eq. 24. 

In the context of mechanics of fracture, it is possible to define 
the so-called stress intensity factors for the delamination in a manner 

H 

analogous to that given in Refs. [4*6] by considering the interlaminar 
stresses ahead of the crack tip along the interface, i.e., 
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(29a) 


Kj - 11m l /ZtT x~^a 2 , (x,0;5 ), 
x-K)^ j-X J J 


U ^ * 

K tt - Urn £ /2? x io,.(x,0|U, 
4,4 xhkO* j«l 0J J 

i 

K T - 11m l fh t x^o,, (x,0;6.) 

- iA X40+ j«l ^ J 


(29b) 


(29c) 


where the superscript a Is omitted, because tractipns, c^, ^4 an<J ° 6 » 
are continuous across the interface. 

While the stress Intensity factors Kj, and describe the 

details of the delaminatlon crack-tip field, the strain energy release 
rate G is also of significant Interest, since this is a quantity 
physically measurable in experiments and mathematically well defined. 
The fracture energy release rate in a delaminated composite may be 
evaluated by using Irwin's virtual crack extension expression [15], 

G " °i + G n + G in 


6 0 

f (a (r,0)[v (k) (68-r,TT) - v (k+1) («P-r,-ir)] 
J 0 y 


■ lim 
60+0 


250 


+ - u (k+1) (SS-r,-ir)] 

+ T (r,0)[w <k) (6S-r,Tt) - w (k+1) C«B-r,-n)])dr , 

yz 


(30) 


where polar coordinates (r,<J>) are used for the convenience of computa- 
tion. The interlaminar stresses, 0 , t , and x in Eq 30, may be 

7 xy yz .. 

obtained from the crack-tip stress field equations such as Eq 27. The 
corresponding displacements are also those of the crack-tip field 
equations obtained in the previous section. 
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5. NUMERICAL EXAMPLES AND DISCUSSION 


The formulation and analysis for the problem outlined in previous 
sections have been programmed into a solution schema suitable for 
numerical computation, For the purpose of illustrating the fundamental 
behavior of the delamination fracture in composite laminates * graphite- 
epoxy systems with symmetric (Q/-9/-0/0) fiber orientation containing 
edge delamination cracks along the Q and -0 ply interface are studied, 
The particular material system and ply orientations are selected here 
because they have been previously investigated in some detail. 

The composite laminate is subjected to a uniform axial extension and 
has a geometry shown in Fi^. 1 with a width-to- thickness ratio 2b/2W 
and uniform ply thickness h^, Delamination cracks of length $ arc 
assumed to emanate from the edges of the composite. Lamina properties 
typical of high-modulus unidirectional graphite-epoxy composite for 
aircraft construction are used in the computation (Table 1). For com- 
posite laminates with the aforementioned laminate geometry and ply 
orientations, several geometric and material symmetry conditions may 
be introduced to simplify the formulation further, The problem, there- 
fore, can be solved very conveniently and accurately. 

5.1 Symmetry and Boundary Conditions, and Further Simplifications 
The symmetric ply orientations and geometry of the composite 
laminate (Fig, 3) lead to the following conditions for displacements: 


on x « b-a, 


(31a) 


3u 3w 3v A 

•Sit m ~— * m — — « ■ 0 

ay 3y 3x 

-¥ 7 “ § 7" 0 ™^-- h 2 - < 31b > 

where the origin of the coordinates is moved to the left tip of the 
delamination* The traction-free boundary conditions on edges and lateral 
surfaces of the composite laminate may be written as 


■if * 0 

xy xz 

on x * -a, 

(31c) 

* t «0 

xy yz 

on y * h^» 

Old) 


Thus, only a quarter of the laminate cross section needs to be considered. 
The boundary conditions, Eqs 31a-d, contain arbitrariness of rigid body 
displacements, which is a characteristic of traction boundary value 
problems . 

Since the eigen solutions, cfj^, satisfy Eq 19 and interface con- 
tinuity conditions Eq 20, we require that a j p ^ satisfy these conditions 

(a) 

also, To determine d£ ' uniquely for Eqs 17 and 18, we further require 
the particular solutions satisfy the following conditions: 


u (p) « 0 

on x " b-cr, 

(32a) 

v (p) « 0 

on y * -h 2 , 

(32b) 

w (p) - 0 

at (0,0,0) . 

(32c) 


Substituting Eqs 17 and 18 into Eqs 19, 20 and 32 gives 

u<“> -I [dk 1) P^ 1) +?£ 1> P^Kb - «), (0-1,2) (33a) 

k*l 
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u 







origiw/'l rw-s ss 
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( 33 b) 

„<«> - 0 , 

( 33 c) 

® * 

( 33 d) 

W 3 a) " 4 ld ^ )+ ^ a) ^ >1, 

( 33 e) 

and 


j i [d' a) + d<“>] .-«£> 

( 34 a) 


( 34 b) 

l [d<“ ) J 0> + d< a) iT<“>) - t (0) , 

k ^ 1 L k H k k M k J xyo’ 

( 34 c) 

Z - Z (df>P< 2 > + df 

k -1 * K K K k -1 K K K K 

( 34 d) 

Z t<»] - f (d« 2) f > + af> I< 2 >], 

k -1 K K k -1 K K K 

( 34 e) 

• j t^&k 0 + Fk ,> «k ,> > + 3 k’ ) fik‘ > + ^ 0> 5 f>>] - 0 . 

( 34 f ) 


Equations 34 a- f give ten linear algebraic equations for twelve real 
unknowns for and u^ . Hence, we may set 

lm[d^ a) J - 0 ( 35 a) 

to reduce the additional degrees of freedom. For symmetric angle-ply 
laminates, it can be shown that Eq 34 f is satisfied identically. 


I 


j! 

!i 

1 


I 


I! 

n 

\\ 


5j 

is 
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Therefore, Instead of using Eq 35a, It may be required 



(35b) 


Since the complete solutions for stresses and displacements must 
satisfy the symmetry and remote boundary conditions, Eqs 31a-d , the 
following relations can be established Immediately to evaluate cr| ' and 



(h) + „(P> . o 
1 1 

<1-1, 5, 6) 

on x * -a, 

(36a) 

(h) (j(p) q 

1 i 

(1-2, 4, 6) 

on y ■ h^ 

(36b) 


and 


3u (h) , 8u (p) . 3v (h >3 V (p > n 3w< h) . aw<P> _ , 

— iy — ■ °‘ W~^W~ " u - 


on x * b~a, 


3v^ . 3v ( ' P ' ) 

— 1 < 4 * ■ — ' 

3x 3x 


3 U < h > 3j/P> „ 3wi h > 3v<P> 

* 3x + 3x “ 3y 


0 , 


(36c) 

on y m -hg* 

(36d) 


By using the eigenfunctions derived previously and the boundary collocation 
metnod, the boundary conditions given in Eqs 36a-d can be matched con- 
veniently in v,he least-square sense. Thus, the eigen solutions for 
cr.jp 1 ) atu j tt (h) can k e determined, Numerical solutions for the problem 
by using the collocation procedure are related to the truncation of 
eigenfunctions and number of collocation stations. Due to space 
limitation, the detailed discussion of solution convergence and accura- 
cy is reported elsehwere [16] • The results presented in this section 

are from collocation calculation, which has a maximum mismatch within 

. o 

one percent deviation from prescribed boundary conditions. A 
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5,2 Strang Singularity for Delamination In Composites 

Now consider a del&mlnatlon lying between 0 and ~0 plies (Fig, 1) 
in a graphite-epoxy composite with ply properties given in table 1, the 
interface continuity and traction boundary conditions along crack sur- 
faces lead to a standard eigenvalue problem for the homogeneous solution, 

ij 

as discussed In Section 3.3, the eigenvalues <5 m obtained from the 
transcendental equation provide basic structures of near-field stress 
and displacement solutions for the delaminafeion problem, the order of 
stress singularity and the asymptotic nature of the crack tip stresses 
depend on the values of 6 m , which satisfy the constraint condition of 
Eq 24, Thus, the eigenvalues corresponding to this restriction are of 
fundamental importance in understanding the delamination failure be- 
havior, For edge delaminated (40) graphite-epoxy composites with 
various fiber orientations, the eigenvalues 6 m , which satisfy the 
aforementioned constraint condition, are found by the present eigen 
analysis and given in Table 2. The stress singularities for an inter- 
face crack between two highly anisotropic laminae are observed to con- 
tain a pair of complex conjugates, 6^ g * -0.541Y, and a constant, 

6 ^ ■ -0,5, This situation is unique and different from that of an 
interface crack between two isotropic or orthotropic media in the sense 
that 5^, 6 2 and 6^ exist simultaneously in the present delamination 
problem, In the degenerated cases such as ±0 ■ 0° and 90°, the composite 
laminates become unidirectional. The delaroination is located in an 
ortho tropic material; the classical inverse square-root singularity 
for crack- tip stresses Is recovered fully. It is noted that the 
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present physical modal and the eigenfunction analysla laid to an oscilla- 
tory stress singularity , as ara tha cases of Interface cracks In iso** 
tropic or orthotropic materials. 

* 

5.3 Asymptotic Stress Field Around Dalaaination 

Complete solutions for delaraination cracks In finite dimensional 
composite laminates are obtainable by using the present Lekhnitskil’s 
complex stress potential formulation and eigenfunction expansion. With 
the aid of the boundary collocation method , the asymptotic stress field 
around a delamination may be expressed in a general form as 


3 

T f /n 

ksr'^ 


7 - 0 . 5+iy 
“k 


-0.5+iY 


n 


'J (k+3) 


4 . (v 

T \M 


jk 


7 -Q.5-lY 

“k 


+ E j (k+3 ) 


_-0.5-iy 
2 > + (F jk 


-0.5 

k 


+ F j (k+3) 



)> 


(j“l» 2 » ..... 16 ) , 


(37) 


where D^, E^, and F^ are known quantities satisfying the following 
relations: 


" E J (k+3 ) ’ D j (k+3) " E ]k’ F Jk " F j (k+3) ’ 


(38) 


to insure being real. More concisely, can be written as 
- L 

Oj * r 2 [A^cos(y An r) + BjSin(y An r) + C^]. * (39) 

For illustrative purposes, the structures of near-field stresses and 
displacements ahead of a delamination (r, 0 ) are given for a 
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<45*/— 45*/-45*/45*) graphite-epoxy laminate with h^ * h 2 « 1 in., 
« * 1 in. and 2b/h ■ 4 subjected to e * * as follows i 

7 * 

a y <r ,0) - [0.04339cos(0. 03434 An r)+Q.39498ain(0, 03434 An r)J 

+ 0 ( 1 ), 

x (r,0) - [0,45347 cos(0. 03434 An r)-0.0498lsin(0. 03434 An r)l 

+ 0 «), 


t xy (ti0) --0. 002449 r 0-5 + 0(1) , 
and 

u (1) (50-r,Tr) - (-0.29576cos[0. 03434 An(60-r)J 

+ 0.01208sin [0.03434 An0$3-r)}} (60~r) 0,5 

- 0.003403(6(3-r) 0,5 + 0(1) , 
v (1) (68-r,ir) « (0.02888cos[0. 03434 An(<50-r)] 

+ 0.70682sin[0. 03434 A n ($0-r)l} (<S0«r) 0,5 

+ 0 ( 1 ), 

w (1 ^(60-r,Tr) - (0.61564 cos[0. 03434 An(60-r)J 

- 0.025155 sin[0. 03434 An(50~r)] } (<$0-r) 0,5 
+ O.OOl6O2(50-r) 0 ” 5 + 0<1), 

where the components of stress are scaled by 10 s * £ z » and the 
displacements by e^. It is noted that the elastic stresses near 


J.-0.5 

(40a) ' 

-0.5 

r 

(40b) 

(40c) 

(4la) 

(41b) 

(41c) 

the 
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delaraination crack tip In a composite laminate possesses the wall known 
oscillatory behavior, and tha displacement field also exhibits an 
oscillatory nature with crack surfaces overlapping each other. As 
first pointed out by Malyshev , et al» (7] and later by England (5] and 
Brdogan [8] for Interfaclal cracks between dissimilar isotropic media, 
the phenomenon of crack surface overlapping is confined to an extremely 
small region and the Interpenetration is not of significance in prac- 
tical terms of fracture mechanics. However, for certain combinations 
of material properties, ply orientations and loading conditions in 
cwfposite laminates, the crack surface contact region has been found to be 
extremely large tl6l . Thus the current model needs to be modified to 
account for the crack surface closure and contact stresses [16]. Studies 
on interface crack closure in dissimilar isotropic media were reported 
recently by Comnlnou [173, Atkinson [18] and Achenback [11] . 


5.4 Delamination Crack Tip Stress Intensity Factors 

Since the Irwin fracture criterion is local in nature and requires 
precise knowledge of the local conditions at the delamination crack tip, 
the stress intensity solutions are obviously of great significance. 
According to the present fracture mechanics theory of composite delaml- 
natlon, stress intensity factors, Kj, and may be evaluated by 

the rigorous analysis described in Section 4* The K^, and 
lead to detailed information of the stress and displacement fields in 
the neighborhood of the delamination crack tip, and may relate to the 
onset of delamination extension upon reaching a critical level. The 
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magnitudes of shown in the formulation depend on the delamlnatlon 

length, ply orientations, laminate geometry, and loading condition?. 

\ 

Consider the (0/-0/-0/0) graphite-epoxy composites with varlou^ 


fiber orientations subjected to uniform axial extension For illus- 
trative purposes, we choose a composite with a width-to-thickness ratio 
2b/ (2h) equal to 8, ply thickness h^ " * 1 in., and delamlnatlon 

length a - 1 in. The mixed mode Kj, and are determined and 


given in Table 3 for various ti'a. It is observed that even though the 
composite laminate is under the simplest loading condition, the delamina- 
tion crack tip response is very complicated due to the complex inter- 
laminar stress distribution, the nonhomogeneity of the solid, the aniso- 
tropic ply properties, and the unusual delamination configuration with 
respect to the loading direction. The out-of -plane tearing mode stress 
intensity factor K.^ caused by interlaminar shear x is about one or 
two orders of magnitude higher than and in general in the lami- 
nates studied. The opening mode stress intensity Kj is also very sig- 
nificant due to the interlaminar normal stresses o . The simultaneous 
presence of K^., and in the delamination problem is unique to 

angle-ply fiber composites, and is not observed in fracture problems 
for bonded dissimilar isotropic media in general. The delamlnatlon be- 
havior is inherently three dimensional in nature; for composites with 
more general laminations, crack geometry and loading conditions, fully 
three-dimensional stress and fracture analyses are essential for ob- 
taining complete information. 

The influence of laminate geometric variables on the delamlnatlon 
behavior is best illustrated by examining the changes of K^., and 






K III tlie ral * tive thickness of upper and lower plies h^/h^ in a 
(45®/— 45*/— 45*/45*) graphite-epoxy composite (with h^ + hj *• W « 2 in.), 
Given the crack length , (fl ■ 1 in.), laminate dimensions (2b ■ 4 in*), 
and the loading condition aa previously, the delamination stress inten- 
sities for various h^/l^'a are shown in Fig* 4. The crack tip tear- 


ing and opening stresses and have a maximum intensification as 


the ply thicknesses h^ become identical, l.e. , h^/hg * 1* The 


however, reaches a minimum due to the reduction of t * It should be 

xy 


noted here that the K^, and depend on material constants of 

all plies as well as the overall geometry* Therefore, the dependency 
of on ply properties is not a simple matter of Identifying them with 
geometric variables, and they may not have the simple physical inter- 
pretation as in the homogeneous case. 


5.5 Strain Energy Release Rates for Delamination 

The equilibrium and stability of delamination are commonly examined 
from an energy rate point of view. The strain energy release rate G de- 
fined in Eq 30 is a quantity characterizing the driving force for delam- 
ination extension. The delamination-growth driving force can be easily 
determined after the establishment of the local asymptotic stress and 
displacement fields. For the edge delamination problem in graphite-epoxy 
composites considered here, the G value may be obtained in a general 
form as 


G — G«, *i* G *4* G 

I II III 


• lira 
60+0 


250 




{A, cos [y in 1 + A 2 sln[y &n ( - ~ - S - ) ] + A 3 > 


(SjLzJL) 

v r ' 


0.5 


dr 


(42) 
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where y ■ Imt^] , Equation 42 has the form similar to the one derived 
previously for an elastic-half space problem by Willis [19], The singu- 
lar integration may be carried out by defining an anlytical function 
with the cut 0 $ x $ 60 


Hz) 



0 


,5+iy 


so that we have 

- r^ ^coafv tn( ^--»—) ]dr - ir«0/(a" Y1T + e YU ), 

Jo r 

‘ 60 ( iL^X ) ° ,5 aln[Y , n( M^r )]dr . 2 yrSB/ (e** Y1T + e Y1T ). 


(43a) 


(43b) 


(43c) 


Thus, the total energy release ratQ G can be determined immediately by 
substituting Eq 43 into Eq 42, Table 4 shows the change of G values with 
ply orientations for the (B/-0/-0/0) graphite-epoxy composites with the 
material properties in Table 1. 

To study the basic nature of delamination extension in angle-ply 
composites, strain energy release rates in the (45 0 /-45 <> /-4S 0 /45 d ) 
graphite-epoxy with various crack lengths are examined. Effects of 
laminate width on delamination crack extension is also investigated. 

The change of total strain energy release rate G with delamination length 
a is given in Tig, 5 to illustrate fundamental characteristics of the 
delamination fracture, For the composite laminates with various 2b/2h’s, 
the G is observed to change with delamination length in a unique manner. 
The maximum energy release rate or crack extension driving force occurs 
at a delamination length approximately equal to one or two ply thick- 
nesses in the composite studied, depending on the (2b/2h) ratio, As 
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Che delamination exceeds this characteristic dimension, G decreases 
monotonically. 

On the baela of fracture mechanics, several important features 

regarding delaminatfon fracture are revealed from the Figure. Assuming 

that the material resistance to delamination growth remains constant 

(i.e., the failure criterion, G - constant, is used), we can iftciedl- 

c 

ately conclude that there exists a critical delamination length asso- 
ciated with the maximum G (for example, a* a 2h for the case b/h m 8) 
for each composite laminate; the word '’critical" means the one that 
experiences stable crack extension at the lowest load. It also indi- 
cates that any Interlaminar edge flaw a Q inherently in the composite, 
which is less than a*, will experience rapidly unstable growth as the 
load or G reaches a critical level, and is anticipated to be arrested 
at a later stage. Any initial delamination greater than a* will ex- 
perience a stable growth under monotonically rising loads; that is, 
there exists an inherently built-in crack arrest mechanism for edge 
delamination. These phenomena predicted by the G - a curve have been 
noted by several researchers conducting experimental and analytical 
studies on the delaraination fracture. The a* may be an important 
quantity in the life prediction for delaminated composite materials 
and structures subjected to static and cyclic loading. 

The delamination strain energy release rate is also a function of 
other geometric variables. For example, G is significantly affected 
by the relative ply thickness h^/hg. In a (45 0 /-45°/-45°/45 0 ) graphite* 
epoxy with a geometry given before, the change of G with h^/hg is given 
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in Fig, 6, whore maximum driving force occurs at h^ * h 2 indicating the 
criticality of the relative ply thickness to delamination fracture in 
composites. 

It is noted here tfyat, oven though the near-field stresses possess 
an oscillatory singularity and may not have the usual significance 
attached to them as in the cohesive (homogeneous) case, the energy re- 
lease rate G is well defined mathematically and physically, and should 
be the quantity of major Interest. The G and its components G^, and 
Gjjj. can be evaluated theoretically and experimentally to provide a 
basic measure of the delamination fracture. 
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6. SUMMARY AMD CONCLUSIONS 


An analytical method for studying delamination is presented in this 
paper. Fundamental nature of edge delamination in advanced fiber composite 
laminates Is examined. Based on the theory of anisotropic elasticity, 
the composite delamination problem is formulated by using Lekhnitskii's 
complex-variable stress potentials and an eigenfunction expansion method. 
Exact orders of the three-dimensional stress singularity in a delamina- 
fclon crack tip region are determined from the eigen analysis. With the 
aid of a boundary collocation technique, complete stress and displacement 
fields in a finite-dimensional, delaminated composite are fully determined. 
Fracture mechanics parameters such as the mixed-mode stress Intensity 
factors and associated energy release rates for edge delamination are 
calculated explicitly. Solutions are obtained for edge-delaminated 
(9/-0/-9/9) angle-ply composites under uniform axial extension, Effects 
of delamination lengths, fiber orientations, lamination and geometric 
variables are studied in detail. 

Based on the information given in the previous sections, the follow- 
ing conclusions may be drawn: 

1, An analytical method based on the theory of anisotropic 
elasticity is successfully developed to study edge de- 
lamlnatlon in angle-ply composite laminates. Formulation 
of the problem is carried out by using Lekhnitskii's complex 
stress functions. Stress singularities for delamination between 
highly anisotropic laminae are obtainable by using an eigen- 
function analysis. The order of delamination crack-tip 
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stress singularity is different from that of an interface 
crack between dissimilar isotropic or orthotropic media 
by the simultaneous presence of three characteristic 
eigenvalues of -0.5+iY, -0,5-iY, and -0.5. 

3. Tue fracture mechanics concept may be extended to delamina- 
tion problems in anisotropic composite laminates by properly 
defining the interlaminar crack-tip stress intensity factors 
such as Eqs 29a-c and strain energy release rates r For 
angle-ply composite laminates, Kj., and always 

occur simultaneously for an edge delamination with being 
one or two orders of magnitude higher than the other two. 

A. Complete stress and displacement fields in a delaminated 
composite may be accurately determined by a combined eigen- 
function expansion and a boundary collocation method, The 
asymptotic solutions are characterized by K ± (1-1,11,111) 
and possess the well known oscillatory behavior. The crack 
surface overlapping could be very large in some composite 
systems with certain combinations of fiber orientations, ply 
stacking sequences, and loading conditions; modifications 
[16] of the current model to include crack closure may be 
needed for these cases. 

5. The crack extension driving force or strain energy release 
rate for edge delamination in composite laminates can be ac- 
curately determined by using Irwin’s crack extension 
concept. Delamination stability in composite laminates 
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under monotonically rising loads can be assessed for any 
inherent interlaminar flaw relative to u.e critical de- 
lamination size a* obtained in the current G - a curve. 

l 


1 
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TABLE 1 


Material’s Constanta for Graphite/ Epoxy Composite Lamina 

E l - 20,0 x 10 6 psl 
E t - 2,1 x 10 6 pal 

°LT “ G Lz " G Tz " °* 85 x 106 P 8 * 

V LT ‘ \z m V Tz - 0>21 
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TABLE 2 


Dominant Stress Singularities for Delamination in (0/-0/-0/O) 

Graphite-Epoxy Composites 


0 

6 1 .2 

ii 

o 9 

- 0.5 

- 0.5 

15° 

- 0,5 ± 0,006421 

- 0,5 

30° 

- 0.5 ± 0.023991 

- 0.5 

45° 

- 0.5 ± 0.034341 

- 0,5 

60° 

- 0.5 ± 0.029421 

- 0.5 

75° 

- 0.5 t 0.015791 

- 0.5 

90° 

- 0.5 

- 0.5 
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TABLE 3 


Stress Intensity Factors K.j* for Edge Delamination In (0/-Q/-O/0) 
Angle-Ply Graphite-Epoxy Composite* Subjected to Uniform Axial Strain 


10 

h 

15° 

0,08645 

30° 

0.2330 

45° 

0.1347 

60° 

0.02025 

75* 

0.006268 


K ++ 

v ++ 

hu 

0,01214 

-4.5588 

0.03330 

-3.6604 

0.01380 

-1,2968 

0.001360 

-0,1775 

0.OO0Z948 

0.0818 


*K, (psl ** /in. ) are scaled by 10 6 e ( 


a m « hj * 1 in. , b 


8 in. 


**For the delamination crack in the first quadrant 


TABLE 4 


Energy Release Rate G fo.r Edge Delamlnatlon in (0/-9/-0/O) Angle -Ply 
Graphite-'Epoxy Composite'}* Under Uniaxial Extension e z 

±0 G/10 6 (psi-in.) 

s 

8.1076 
4.0506 
0.5740 
0.0138 
0.0036 

■ hg 1 in., b ■ 8 in. 


30° 

45’ 

60° 

75° 
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FREE EDGE INTERFACE 



DELAMINATION GEOMETRY AND LAMINATE CONFIGURATION 


ORIGINAL PAGE !S 

OF POOR QUALITY 






FIG. 4 STRESS 
( 45 % 45 * 
RATIO I 




STRAIN ENERGY RELEASE RATE G vs DELAMINATION CRACK LENGTH a IN (45%-45 0 /-45 0 /45°) GRAPHITE-EPOXY 
COMPOSITES WITH VARIOUS LAMINATE WIDTHS b ’s (hi = h 2 = I ia.) 
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FIG. 6 ENERGY RELEASE RATE G FOR EDGE DELAMWATION CRACK IN (4 5% 45^-45^4^} 


APPENDIX I 


Materials Parameters for Tielamlnatlon Crack-Tip Eigenvalues 
in Cross-Ply Composite Laminate 


<$ * n 

n 


or 


* (n - ^) ± 7 ^ An{[b + (b 2 - 4a 2 ) 2 ]/(2a)> 


n 

where 


a - + M 3 (k) M^ k+1 ^ + M^ k) M^ k+1 ^ 


+ M^ k) + M^ k+1 ) 
4 4 


b » - 2M.[ k ^ k+1 ^ - 2M<j k) M 

- 2 M^ k+1 ’ ) - 2Mp^ M^ k+1 ^ - 2M, ( ' k+1 ^Mj k+1 ' ) 
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